1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 4000000 iteration;
- 3000000 burn-in;
- 1000 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.00
## LDevsum 1.02 1.08
## dh0 1.00 1.02
## dh1 1.00 1.01
## dl0 1.05 1.23
## dl1 1.02 1.07
## dl2 1.05 1.20
## dl3 1.00 1.00
##
## Multivariate psrf
##
## 1.04
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1211.71099 | 22164.3356 |
| DIC3 | 1149.27410 | 22188.2202 |
| PWAIC | 41.02646 | 242.9939 |
| WAIC | 1175.11970 | 22209.9391 |
2. Separate Modeling of High-Level
2.1 Trace plots for convergence check
The current MCMC setting is:
- 4000000 iteration;
- 3000000 burn-in;
- 1000 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.000 1.01
## dh0 0.999 1.00
## dh1 0.999 1.00
##
## Multivariate psrf
##
## 1
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| H0a | |
|---|---|
| DIC | 1208.05552 |
| DIC3 | 1149.19337 |
| PWAIC | 40.82378 |
| WAIC | 1174.68288 |
3. Separate Modeling for Low-level
3.1 Trace plots for convergence check
The current MCMC setting is:
- 4000000 iteration;
- 3000000 burn-in;
- 1000 thinning.
3.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## LDevsum 1.020 1.095
## dl0 1.063 1.223
## dl1 1.009 1.044
## dl2 1.058 1.198
## dl3 0.999 0.999
##
## Multivariate psrf
##
## 1.03
3.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
3.4 WAIC results
| L4 | |
|---|---|
| DIC | 22161.3942 |
| DIC3 | 22186.6573 |
| PWAIC | 244.2203 |
| WAIC | 22207.8939 |